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PLASTIC  STRESS-STRAIN  RELATIONS  BASED  ON 
LINEAR  LOADING  FUNCTIONS^ 

o 

J.  L.  Sanders,  Zr, 


ABSTRACT 


The  stress-strain  relations  of  Incremental 
theories  of  plasticity  have  recently  been  general- 
ized by  Koiter  vho  introduced  the  use  of  any  number 
of  loading  functions  instead  of  only  one.  The 
present  paper  is  concerned  with  incremental  theory 
based  on  any  number  of  linear  loading  functions. 

In  this  case  the  resulting  plastic  strain  incre- 
ments are  integrable  in  a restricted  sense  and  the 
stress-strain  relations  can  be  given  in  a form 
partially  resembling  deformation  theories.  The 
slip  theory  of  Batdorf  and  Budiansky  falls  within 
the  class  of  theories  which  are  discussed  in  this 
paper. 


^he  results  presented  in  this  paper  were  obtained  in  the  course 
of  research  sponsored  by  the  Office  of  Naval  Research  under 
Contract  N7onr-35C01, 
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Research  Associate  of  Applied  Mathematics,  Brown  University. 
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IntroduQtlon.  Within  the  last  few  years  the  development  of  plastic 
stress -strain  relations  for  strain  hardening  materials  has  received 
increasing  attention  from  workers  in  the  field  of  applied  mechanics 
[1]|[2],[3],[4],^  a knowledge  of  the  stress-strain  curve 
for  a material  in  simple  tension  or  compression  is  not  in  itself 
sufficient  for  predicting  the  behavior  of  a material  under  more 
complex  types  of  loading.  At  least  two  problems  iinmedlatly  pre- 
sent themselves  when  one  considers  loading  which  may  Involve  all 
six  components  of  stress.  These  are  the  generalization  of  the  idea 
of  a yield  point  to  the  idea  of  a yield  surface  in  a stress  space 
of  six  dimensions  and  the  formulation  of  a plastic  stress-strain 
relation  for  polyaxial  loading. 

During  the  history  of  the  subject  these  problems  have 
been  attacked  in  two  more  or  less  distinct  ways.  One  method  of 
approach  is  to  try  to  deduce  the  behavior  of  a bulk  piece  of  metal 
from  the  behavior  of  the  individual  crystals  making  up  the  aggreg- 
ate, A recent  example  of  results  obtained  by  pursuing  this  method 
is  the  slip  theory  of  plasticity  [ 2 ] , The  other  method  of  approach 
is  to  try  to  deduce  from  a few  reasonable  assumptions  what  the 
general  form  of  the  yield  surface  and  the  stress-strain  relation 
must  be  without  inquiring  very  deeply  into  the  structure  of  metals 
but  rather  proceeding  in  a phenomenological  manner  [ 3 ] • 

Results  obtained  by  one  method  often  suggest  research 
to  be  done  by  the  other  method  which  tends  to  unify  the  whole  sub- 
ject. The  earliest  known  yield  condition,  the  Tresca  maximum  shear 

^Numbers  in  brackets  refer  to  the  bibliography  at  the  end  of  the 
paper. 
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condition  [ 5 } | when  expressed  in  geometrical  form  as  a yield  sur- 
face in  stress  space  possesses  "corners".  The  initial  yield  sur- 
face according  to  slip  theory  is  the  Tresca  yield  surface  which 
arises  as  an  envelope  of  planes  in  stress  space.  The  yield  surface 
after  work  hardening  according  to  slip  theory  is  also  an  envelope 
of  planes  and  naturally  possesses  corners.  Drucker  [3]  allows  for 
corners  in  his  general  discussion  of  yield  surfaces  from  a phenom- 
enological point  of  view.  The  relation  between  corners  in  the  yield 
surface  and  linearity  of  the  stress-strain  relation  was  investigated 
experimentally  by  Stockton  [ 6 ] . The  yield  surface  is  closely 
connected  with  a function  called  the  loading  function  which  appears 
in  the  stress-strain  relation.  The  geometrical  expression  of  the 
loading  function,  the  loading  surface,  is  identical  to  the  yield 
surface  if  there  is  only  one  loading  function  in  the  stress-strain 
relation;  if  the  stress-strain  relation  is  based  on  many  loading 
functions  the  yield  surface  is  the  boundary  of  those  points  in 
stress  space  which  represent  elastic  behavior  with  respect  to  all 
the  loading  surfaces.  Obviously  corners  in  the  yield  surface  may 
arise  even  though  all  the  loading  surfaces  are  smooth.  Stress- 
strain  relations  based  on  many  loading  functions  were  suggested  by 
Koiter  who  also  investigated  some  of  their  properties  [ 4 ] , Al- 
though slip  theory  was  derived  through  physical  reasoning  Koiter 
states  that  it  may  be  regarded  as  a theory  of  plasticity  based  on 
an  infinite  number  of  plane  loading  surfaces. 

An  interesting  feature  of  slip  theory  is  that  the 
incremental  form  of  its  stress-strain  relation  is  Integrablo  in  a 
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restricted  sense  which  resembles  deformation  theories  In  this  re- 
spect and  Is  unlike  most  Incremental  theories.  The  present  paper 
Is  concerned  with  the  use  of  many  linear  loading  functions  In  con- 
structing stress-strain  relations.  It  Is  shown  that  in  some  cases 
the  resulting  incremental  stress-strain  relation  is  Integrable  in 
a restricted  sense.  Means,  suggested  by  the  behavior  of  a single 
crystal,  of  including  cross  effect  and  Bauschingor  effect  in  the 
stress-strain  relation  are  also  discussed. 

A Single  Plane  Loading  Surface.  In  the  formulation  of  a plastic 
stress-strain  relation  it  is  commonly  assumed  that  the  strain  may 
be  separated  Into  an  elastic  and  a plastic  part  which  arc  additive. 
This  assumption  is  made  in  the  following  treatment  and  the  plastic 
strain  only  is  the  subject  of  discussion.  It  is  further  assumed  in 
the  following  that  the  plastic  part  of  the  strain  may  be  separated 
into  many  parts  each  associated  with  a loading  function. 

It  Is  convenient  to  phrase  the  discussion  in  geo- 
metrical terms  and  spcalf:  of  loading  surfaces,  loading  paths,  strain 
increment  vectors  and  so  forth.  Imagine  a euclidean  space  with 
rectangular  coordinates,  one  coordinate  for  each  component  of 
stress.  Each  stress  state  is  represented  by  a point  in  this  space. 
For  convenience  wo  also  represent  strain  in  the  same  space  with 
each  component  of  the  strain  acting  as  coordinate  on  the  same  axis 
as  the  corresponding^  component  of  stress.  The  loading  history  of 
a piece  of  strain-hardening  material  which  has  been  subjected  to  a 
varying  homogeneous  state  of  stress  is  representable  in  this  space 
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as  a loading  path.  Wc  now  want  to  consider  tho  contribution  of  a 
single  piano  loading  surface  to  tho  plastic  strain. 

The  stress-strain  relation  to  bo  described  in  the 
following  is  consistent  in  that  it  satisfies  all  the  mathematical 
conditions  which  have  been  discovered  by  phenomenological  studios 
to  apply  to  incremental  stress  theories  of  plasticity.  The  proof 
that  in  tho  particular  case  of  a plane  loading  surface  the  in- 
cremental stress -strain  relation  is  intograble  is  given  in  tho 
appendix,  Tho  stress-strain  relation  in  its  integrated  form  is 
presented  here. 

Fig,  1 shows  a symbolic  representation  of  the  loading 
plane  and  the  loading  paths  of  two  hypothetical  experiments.  As 
loading  proceeds  from  0 to  P the  material  behaves  elastically.  At 
P plastic  straining  begins  and  continues  as  long  as  the  loading  path 
is  directed  outward  from  the  loading  plane.  Thus  during  the  part 
Q'R  of  tho  loading  path  OPQ'R  no  plastic  straining  occurs.  One  may 
think  of  the  loading  plane  being  pushed  out  parallel  to  itself  as 
loading  proceeds.  The  total  plastic  strain  suffered  by  this  hypo- 
thetical material  depends  only  on  the  distance  the  loading  plane 
has  been  pushed  out.  Thus  the  total  plastic  strain  is  the  same  for 
either  loading  path  OPQ  or  OPQ'R  and  is  a function  of  d.  Tho  direc- 
tion of  the  plastic  strain  vector  is  perpendicular  to  tho  loading 
plane. 

Two  Plane  Loading  Surfaces. 

Several  novel  features  of  tho  use  of  many  plane  load- 
ing surfaces  can  be  illustrated  most  simply  by  examining  the  case  in 
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which  there  arc  only  two  piano  loading  suri'acos  and  f2«  Fig,  2 
is  a schematic  drawing  which  shows  how  stross  space  is  divided  into 
four  regions  by  two  plane  loading  s^irfacos.  If  the  loading  path 
remains  in  region  I no  plastic  straining  occurs.  If  the  loading 
path  only  penetrates  region  II  (or  IV)  then  only  one  of  the  two 
pianos  is  affected  and  we  have  the  case  treated  in  the  last  section. 
If  however  the  stross  path  enters  region  III  then  both  loading  planes 
are  affected  as  shown  in  Fig.  3. 

Suppose  the  loading  haa  reached  a point  such  as  P in 
Fig.  3 and  now  a small  increment  of  stress  do.,  is  added  to  the 

^ J 

load.  Th'^re  are  several  cases  to  consider  depending  on  the  direction 
of  the  stress  Increment  vector.  Pig.  4 shows  the  region  of  stross 
space  in  the  neighborhood  of  P.  There  arc  four  regions  around  P 
into  which  the  stress  increment  vector  may  extend.  If  dc^j  extends 
into  region  1 there  is  no  additional  plastic  strain,  consequently 
region  1 may  be  called  the  elastic  region.  If  do.  extends  into 

* J 

either  region  2 or  4 one  of  the  two  pianos  is  affected  and  the  other 
is  not.  If  do^j  extends  into  region  3 both  loading  planes  are 
affected  so  region  3 may  be  called  the  region  of  total  loading. 

The  strain  increment  is  taken  to  be  the  sum  of  the  contribution 
from  fj^  and  f2  as  though  each  acts  independently. 

The  boundary  of  the  clastic  region  is  the  yield  sur- 
face which  is  shown  in  heavy  lines  in  Pig.  4,  In  the  case  under 
consideration  a "corner"  has  developed  in  the  yield  surface  at  the 
load  point  P.  If  the  load  point  wore  on  any  straight  portion  of 
the  yield  surface  the  strain  increment  vector  would  of  course  be 
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normal  to  tho  yield  surface.  At  a corner  such  as  in  Fig.  4 there 
is  no  normal  to  tho  yield  surface,  Howovor,  the  strain  increment 
vector  dE^j  is  tho  vector  sum  of  tho  contribution  from  both  load- 
ing surfaces.  Therefore  at  a corner  tho  strain  Increment  vector 
has  some  determinate  direction  between  tho  normal  to  and  the 
normal  to  fg  (see  Fig.  5). 

Tho  mathematical  expression  of  tho  stress-strain 

relation  which  we  are  considering  contains  the  differentials  of 

stress  ddj,j  to  the  first  power  only  and  thus  is  linear  in  some 

respects.  If  the  stress  increment  vector  do^^j  causes  a plastic 

strain  de|^  then  ^dc^^j  causes  twice  as  mucli  plastic  strain,  2de^j. 

However  -do  mi^ht  not  cause  any  plastic  strain  at  all  because  if 
ij 

the  vector  do^^  falls  in  the  region  of  total  loading  then 
falls  in  tho  elastic  region,  Hooke’s  law  for  the  elastic  part  of 
the  strain  is  truly  linear  and  does  not  possess  this  peculiarity. 

The  plastic  stress-strain  relation  is  non-linear  because  tho  prin- 
ciple of  superposition  does  not  apply,  A plastic  stress-strain 
relation  which  involves  even  only  one  loading  function  is  non-linear 
in  this  respect. 

The  non-linearity  for  multiple  loading  surfaces  is 


of  a less  trivial  sort.  Suppose  for  instance  that  do 


ij(C) 


do 


iJ(A) 


^‘’iJ(C)  acting  by 


itself  causes  non-zero  plastic  strain  increments  der^. . . , de?.,  » 

ij(A)  ’ ij(3) 

a^ii  rospectivoly.  In  general  it  docs  not  follow  that 

® dfiP  , - + FiS*  6 illustrates  a case  in  which 

ijvc;  ij(A;  ij(B; 
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superposition  does  not  hold.  The  stress  increment  do-.>.»  is  applied 

ij(A) 

followed  by  Clearly  the  resulting  plastic  strain  is  not 

the  same  as  that  produced  by  acting  alone  because  the  load- 

ing plane  is  moved  out  more  in  the  first  Instance  than  in  the 

second.  Even  the  order  in  which  do..,.,  and  d<j..,„v  are  apnlied 

lj(A;  IjQd; 

makes  a difference.  Experimental  evidence  for  the  existence  of 
corners  in  the  yield  surface  has  recently  been  found [ 6 ] , 


Many  Plane  Loading  Surfaces.  There  is  no  difficulty  in  generaliz- 
ing the  previous  results  by  introducing  any  finite  or  infinite 
number  of  plane  loading  surfaces.  If  there  are  an  Infinite  number 
of  loading  surfaces  the  stress-strain  relation  may  involve  an  in- 
tegration, The  set  of  all  those  points  in  stress  space  which  are 
elastic  with  respect  to  all  the  loading  planes  constitutes  the 
elastic  region.  The  boundary  of  the  elastic  region  is  the  yield 
surface.  Figs.  7 and  8 show  symbolically  the  elastic  region  and 
yield  surface  when  the  number  of  plane  loading  surfaces  is  finite 
or  infinite  respectively. 

As  loading  proceeds  the  loading  planes  get  pushed  out 
more  or  less  depending  on  the  loading  path,  but  at  any  instant  of 
loading  there  is  a well  defined  yield  surface.  The  property  of  in- 
tegrability  for  a single  plane  loading  surface  leads  to  the  follow- 
ing simple  rule  in  the  case  of  many  plane  loading  surfaces i 


The  total  plastic  strain  is  the  same 
for  all  loading  paths  which  result  in 
the  same  yield  surface. 
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Pig,  9 illustrates  a ease  in  which  two  different  loading  paths  A 
and  B produce  the  same  total  plastic  strain.  If  tho  loading  had 
proceeded  only  as  far  as  P the  total  strain  for  paths  A and  B 
would  bo  different  but  continuing  the  loading  to  Q wipes  out  this 
difference, 

Tho  gonorallzatlon  of  the  notion  "region  of  total 
loading"  to  tho  ease  of  many  plane  loading  surfaces  can  bo  done  in 
the  following  way.  Suppose  loading  has  procooiod  along  a straight 
line  loading  path  into  the  plastic  range  and  reached  a point  P 
(Fig,  10).  The  yield  surface  is  made  up  of  part  of  the  Initial 
yield  surface  and  of  part  of  the  envelope  of  the  loading  planes 
which  pass  through  P and  are  tangent  to  the  initial  yield  surface. 

The  part  of  this  conical  envelope  which  extends  outward  from  P 
bounds  the  region  of  total  loading.  All  loading  planes  previously 
pushed  out  in  reaching  P are  further  pushed  out  by  a stress  in- 
crement dd  whose  vector  drawn  from  P falls  in  the  region  of  total 
loading. 

The  straight  line  loading  path  (radial  loading)  is  a 
special  case  of  a more  general  loading  path  traced  out  by  a load 
point  P which  always  proceeds  into  the  region  of  total  loading. 

For  such  loading  paths  (and  only  for  such  loading  paths)  the  interest- 
ing result  has  been  obtained  that  J2  deformation  theory  coincides 
with  a particular  incremental  theory  based  on  an  infinite  number 
of  plane  loading  surfaces.  The  details  are  left  to  a later  paper. 

Further  Remarks  on  Plane  Loading  Surfaces.  In  all  the  foregoing 
discussion  it  has  been  assumed  that  the  loading  planes  act  entirely 
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Independently.  Fig,  11  Illustrates  a case  in  which  the  loading 
planes  act  interdependently.  The  loading  plane  has  been  pushed 

out  to  the  position  indicated  by  the  dotted  line.  The  other  planes 
also  move  according  to  some  law  as  a consequence  of  the  motion  of 
f^.  This  type  of  behavior  is  certainly  more  complicated  than  in- 
dependent action  but  the  phenomenon  of  latent  strain  hardening  of 
a single  crystal  suggests  such  behavior.  In  a metal  crystal  there 
is  a certain  preferred  set  of  planes  called  slip  planes  so  oriented 
with  respect  to  the  crystal  lattice  that  slip  is  most  likely  to 
occur  in  them  rather  than  in  any  others.  Within  these  planes  there 
are  several  preferred  directions  of  slip.  Thus  there  is  some 
definite  number  of  slip-plane  slip-direction  combinations  or  slip 
systems  within  a given  crystal.  When  a given  slip  system  has  been 
activated  it  becomes  more  resistant  to  further  slipping,  that  is 
it  becomes  strain  hardened.  However  the  slip  systems  do  not  strain 
harden  independently,  V/hen  any  one  slip  system  is  strain  hardened 
so  are  the  others.  The  condition  that  the  shear  stress  on  a given 
slip  system  remain  constant  is  that  a certain  linear  combination 
of  the  stress  components  remain  constant.  The  locus  of  points  in 
stress  space  satisfying  this  condition  is  of  course  a plane.  Per- 
haps a theory  describing  the  macroscopic  plastic  strains  of  a single 
crystal  could  be  formulated  based  on  a finite  number  of  plane  load- 
ing surfaces  acting  interdependently.  Even  though  it  is  not  at  all 
obvious  how  to  take  the  step  from  a single  crystal  to  a crystal 
aggregate  the  behavior  of  a single  crystal  is  very  suggestive. 

Fig,  12a  shows  a possible  yield  surface  after  an 
experiment  performed  on  a material  whose  initial  yield  surface  is 
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2 

J2  « k . The  material  has  been  stressed  In  simple  tension  in  the 
X direction  beyond  the  elastic  limit.  The  yield  stresses  for 
tension  and  compression  in  the  y direction  have  been  altered 
(cross  effect)  and  the  yield  stress  for  compression  in  the  x direc- 
tion has  also  been  altered . (Bauschinger  effect)  [ 7]  . Fig.  12b 
shows  the  results  of  the  same  experiment  performed  on  a hypothetical 
material  which  obeys  a stress-strain  relation  based  on  an  Infinite 
number  of  plane  loading  surfaces  which  act  independently.  The 

p 

initial  yield  surface  for  this  material  is  also  J2  *=  k , No  cross 
effect  or  Bauschinger  effect  in  the  sense  of  a lowering  of  the 
yield  point  in  compression  followin,  loading  in  tension  is  exhibited 
by  such  a material.  Fig.  12c  shows  the  result  of  the  same  experi- 
ment performed  on  another  hypothetical  material  obeying  the  simple 
3 2 rule.  The  material  of  12c  could  also  be  regarded  as  obeying  a 
stress-strain  relation  based  on  an  infinite  number  of  plane  loading 
surfaces  in  which  the  planes  always  move  in  such  a way  as  to  envel- 
op a surface  3^  = constant.  The  property  of  integrability  does  not 
necessarily  hold  for  stress-strain  relations  based  on  many  plane 
loading  surfaces  which  act  interdependently.  The  properties  of 
such  stress-strain  relations  have  only  begun  to  be  investigated, 
but  the  phenomenon  of  latent  strain  hardening  exhibited  by  individ- 
ual metal  crystals  as  well  as  the  evidence  from  experiments  on 
polyorystalline  materials  suggests  further  investigations, 

Many  plate  buckling  experiments  agree  well  with  a 
theory  in  which  the  hypothetical  material  of  Fig,  12b  is  taken  as 
the  mathematical  abstraction  of  the  real  material.  Unfortunately 
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other  kinds  of  experiments  show  that  this  material  is  not  accept- 
able as  a good  approximation  to  a real  material.  Recent  investiga- 
tions [ 8 ] on  the  stability  of  non-linear  systems  indicate  that  the 
classical  buckling  criterion  is  not  applicable  in  the  plate  buckling 
problem.  Those  investigations  show  that  small  and  unavoidable 
initial  imperfections  can  have  a pronounced  effect  on  the  load 
carrying  capacity  of  a plate  and  yet  not  cause  much  scatter  in  the 
test  data.  If  those  conclusions  are  correct  plate  buckling  experi- 
ments are  not  decisive  in  the  choice  between  various  possible  stress- 
strain  relations. 


! 
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APPENDIX 


Consider  the  following  stress-strain  relatloni 


de^  s 

IJ 


0(f)  df  df  > o 

ss  O O 


(1) 


In  which  f,  the  loading  function,  is  a function  of  stresses  alone. 
Provided  the  surface  f = constant  bo  sufficiently  smooth  and  non- 
convex  toward  the  origin  in  stress  space,  this  is  the  most  general 
stress-strain  relation  involving  only  one  loading  function  which 
belongs  to  the  class  of  stress  theories  of  plasticity.  If  f is 
a linear  combination  of  the  stresses  then  (1)  may  be  integrated  to 
give  the  components  of  total  strain  in  terms  of  the  stresses.  The 
resulting  relation  is  of  the  formi 


P _ 


ij 


= H(f)  a 


iJ 


(2) 


where  dH  » Odf  and  o are  direction  numbers  of  the  normal  to  the 

loading  plane.  This  may  be  verified  by  taking  the  differential  of 

(2)  to  obtain  (1),  Conversely  the  only  intcgrablc  stress-strain 

relations  of  the  form  (1)  are  those  for  which  f = constant  is  a 

plane  surface  in  stress  space.  This  may  be  proved  by  the  following 

considerations.  If  de^j  is  intograblo  then  there  exists  a family 

of  surfaces  in  stress  space  such  that  is  constant  on  each 

J 

surface  of  the  family.  Supnoso  that  a point  moves  on  a surface 

f a constant}  then  according  to  (1)  df  * o and  de^[*  = o so  the 

P 

surface  f = constant  must  also  bo  a surface  eT,  * constant,  hence 

* J 
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Is  a function  of  f alone.  This  implies  that  is  a function 

^ij 

of  f alone,  sayi 


'ij 


(3) 


Now  J£  are  the  direction  numbers  of  the  normal  to  the  surface 

f = constant  and  (3)  implies  that  the  normal  has  a constant  direc- 
tion; hence  the  surface  f = constant  is  a plane.  All  that  has  been 
proved  is  that  f must  be  a function  of  a linear  combination  of  the 
stresses,  however  there  is  no  loss  in  generality  in  taking  f itself 
to  be  linear  in  the  stresses  since  any  additional  arbitrariness  may 
be  absorbed  into  the  factor  G which  occurs  in  (1), 


All-97 


15 


Bibliography 


1.  Handelman,  G.  H«:  Lin,  C.  C.;  Prager.  W* ; "On  the  Mechanical 
Behavior  on  Metals  in  the  Straln-Haraenlng  Range"  Quar.  Appl. 
Math.  V. 4(1947)  pp.  397-407. 

2.  Batdorf,  S.  B, ; Budiansky,  B.;  "A  Mathematical  Theory  of 
Plasticity  Based  on  the  Concept  of  Slip".  N.A.C.A.  TN  1871 
(1949). 

3.  Drucker,  D.  C.;  "A  More  Fundamental  Approach  to  Plastic  Stress- 
Strain  Relations."  Proceedlngsof  the  First  National  Congress 
of  Applied  Mechanics.  ASME  (1951)  pp.  487-491. 

4.  Kolter,  W.  T.;  "Stress-Strain  Relations,  Uniqueness  and  Varia- 
tional Theorems  for  Elastic-Plastic  Materials  with  a Singular 
Yield  Surface."  Quar.  Appl.  Math.  V.  11  (Oct.  1953). 

5.  Tresca,  H.;  "Memoire  sur  I'ecoulement  des  corps  solides  soumis 
a de  fortes  pressions."  Coroptes  rendus  de  1' academic  des 
sciences  Paris.  V.59  (1864)  pp,  754-758. 

6.  Stockton,  F.  D.5  "Experimental  Evidence  of  Non-Linearity  in 
Plastic  Stress-Strain  Relations."  Report  All-88  to  ONR  (1953). 
Graduate  Division  of  Appl.  Math.,  Brown  University, 

7.  Edelman,  F.  and  Drucker,  D.  C.;  "Some  Extensions  of  Elementary 
Plasticity  Theory,"  Jour.  Franklin  Inst,  V.  251  N.6  (Jun,  1951) 

pp.  581-605. 

P,  Onat,  E,  T,  and  Drucker,  D,  C,  : "Inelastic  Instability  and 
Incremental  Theories  of  Plasticity,"  Jour.  Aero,  Sci,  V,  20 
(Mar.  1953)  PP.  I8I-I86. 


FIG.  M 


Tntcrdepf*;v.e-:t  i.lu'.e  ioiding  3urr?.ce.*? 


APPROVED  DISTRIBUTION  LIST  FOR  UNCLASSIFIED  TECHNICAL  REPORTS 

lasiMd  by 
BROWN  UNIVERSITY 
Contract  N7onr-358,  T.O.  X 
NR  OM-1  032 


Office  of  Naval  Research 
Washington  25y  D.  C« 

H-2  Attn:  Mathematics  Branch  (Code  432) 
M-1  Mechanics  Branch  (Code  438) 
M-1  Physics  Branch  (Code  421) 
M-1  Metallurgy  Branch  (Code  423) 


M-2 


M-1 


M-1 


M-1 


M-1 


M-18 


M-9 


M-2 

M-1 

M-1 

M-1 


Conmandlng  Officer 

Office  of  Naval  Research  Branch  Office 
150  Causeway  Street 
Boston)  Massachusetts 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 

346  Broadway 

New  York,  Now  York 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 
844  North  Rush  Street 
Chicago  11,  Illinois 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 

1000  Geary  Street 

San  Francisco  9,  California 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 
1030  East  Green  Street 
Pasadena  1,  California 

Of  fleer- In-Char  g e 
Office  of  Naval  Research 
Navy  #100 
Fleet  Post  Office 
New  York,  New  York 

Director 

Naval  Research  Laboratory 
Washington  20,  D.  C, 

Attn:  Scientific  Information  Division 
Library  (Code  2021) 

Applied  Mathematics  Branch  (Code  3830) 
Shock  and  Vibrations  Section  (Code  3850) 
Structures  Branch  (Code  3860) 


M-2 

R-1 

M-1 


M-2 

M-1 


M-2 


M-1 

R-1 

R-1 


R-2 


M-1 

M-1 


R-1 


M-1 

M-1 


M-1 

M-1 


Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn*  Code  364  (Technical  Library) 

Code  423  (Underwater  Explosion  Research; 
Code  442  (Scientific  Section,  Design) 


David  Taylor  Model  Basin 
Carderock.  Maryland 
Attn*  Library 

Structural  Mechanics  Division 


Naval  Ordnance  Laboratory 
White  Oak.  Silver  Spring  19,  Maryland 
Attn*  Library 

Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.  C* 

Attn*  AER-TD-414 

Materials  Branch 
Design  Elements  Division 

Bxireau  of  Yards  and  Docks 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Director,  Research  Division 

Commander 

Norfolk  Naval  Shipyard 
Norfolk,  Virginia 

Attn:  Technical  Library  (Code  243 A) 
UERD  (Code  290) 


Superintendent 

Aeronautical  Structures  Laboratory 
Building  600,  Naval  Air  Experimental  Station 
Philadelphia  12,  Pennsylvania 

Attn:  Experimental  Structures  Section 

Office,  Assistant  Chief  of  Staff,  G4 
The  Pentagon 
Washington,  D.  C« 

Attn:  Research  and  Development  Division 

The  Chief.  Armed  Forces  Special  Weapons  Project 
Department  of  Defense 
P.  0.  Box  2610 
Washington,  D.  C« 

U.  S.  Army  Arsenal 
Watertown  72,  Massachusetts 
Attn:  Dr.  R.  Beeuwkes 
Mr.  J.  Bluhm 


3 


Prankford  Arsenal 
Pitraan-Dunn  Laboratory 
Philadelphia  37,  Pennsylvania 
M-1  Attn*  Dr.  Herbert  I.  Pusfeld 

Picatinny  Arsenal 
Dover,  New  Jersey 
M-1  Attn:  Drt  L*  Gilman 

Commanding  General 
V/right  Air  Development  Center 
V/right-Patterson  Air  Force  Base 
Dayton,  Ohio 
MJ+  Attm  WCACD 

Department  of  Commerce 
Office  of  Technical  Service 
V/ashington  25,  D«  C, 

M-1  Attn:  Library  Section 

National  Advisory  Committee  for  Aeronautics 
1724  F.  Street  NW 
V/ashington  25*  D.  C« 

M-1  Attn:  Chief  of  Aeronautical  Intelligence 

National  Advisory  Committee  for  Aeronautics 
Langley  Aeronautical  Laboratory 
Langley  Field,  Virginia 
M-1  Attn:  Library 

National  Advisory  Committee  for  Aeronautics 
Lewis  Flight  Propulsion  Laboratory 
Cleveland  Airport 
Cleveland  11,  Ohio 
M-1  Attn:  Library 

National  Bureau  of  Standards 
Washington,  D.  C* 

M-1  Attn:  Dr.  W.  R.  Ramberg 

Director  of  Research 
Sand la  Corporation 
Albuquerque,  New  Mexico 
M-1  Attn:  Dr.  S.  C.  Right 

Brooklyn  Polytechnic  Institute 
85  Livingston  Street 
Brooklyn,  Now,  York 
R-1  Attn*  Dr.  N.  J.  Hoff 

R-1  Dr.  H.  Relssner 

M-1  Dr.  P.  G.  Hodge,  Jr. 

M-1  Dr.  F.  S*  Shaw  (Dept.  Aero.  Engr.  & Appl.  Mech. ) 

Brown  University 
Providence  12,  ^ode  Island 

M-1  Attn:  Chairman,  Graduate  Division  of  Applied  Mathematics 


4 


R-1 

R-1 

R-1 


M-1 

R-1 

R-1 

M-1 


R-1 

R-1 


R-1 

R-1 

M-1 

R-1 

R-1 


M-1 

R-1 


M-1 

R-1 


R-1 


M-1 

M-1 

M-1 

R-1 


R-1 

R-1 


M-1 


California  Institute  of  Technology 
Pasadena,  California 

Attnt  Dr*  J*  G*  Kirkwood 
Dr.  Pol  Duwez 
Dr*  H*  S«  Tslen 

University  of  California 
Berkeley,  California 
Attni  Dr*  J*  E*  Dorn 
Dr*  H«  Hultgren 
Dr*  G«  C*  Evans 
Dr*  C.  F*  Garland 

University  of  California 
Los  Angeles,  California 

Attn:  Dr*  I*  S.  Sokolnlkoff 

Dr*  D.  Rosenthal 

Carnegie  Institute  of  Technology 
Pittsburgh,  Pennsylvania 
Attn:  Dr*  J.  S*  Koehler 

Dr*  G*  H*  Handelraan 
Dr.  E.  Saibel 
Dr*  II*  J.  Greenberg 
Dr.  E.  D'Appolonla 

Case  Institute  of  Technology 
Cleveland,  Ohio 

Attn:  Dr*  U.  M.  Baldwin,  Jr* , Metals  Research 

Dr*  0.  Hoffman 

Catholic  University  of  America 
Washington,  D.  C. 

Attn:  Dr*  F.  A.  Biber stein 

Dr*  K*  Hertzfeld 

University  of  Chicago 
Chicago,  Illinois 

Attn:  Dr*  C.  S.  Barrett 

Columbia  University 
New  York,  New  York 

Attn:  Dr*  R.  D*  Nindlln 

Dr*  H*  Bleich 
Dr*  B*  A*  Bo  ley 
Dr*  M*  Gensamor 

Cornell  University 
Ithaca,  New  York 

Attn:  Dr*  H*  S*  Sack 
Dr*  A Kant r owl tz 

University  of  Florida 
Gainesville,  Florida 

Attn:  Dr*  C.  G*  Smith 


Labora- 

tory 


I 


f 


5 


R-l 

M-1 


R-l 

M-1 

R-l 


M-1 

R-l 

R-l 

R-l 

M-1 


M-1 


R-l 


R-l 

R-l 


M-1 

M-1 


R-l 


R-l 

R-l 

R-l 

R-l 

R-l 

R-l 

R-l 

R-l 

M-1 


Harvard  University 
Cambridge  381  Massachusetts 

Attm  Dr.  P.  Birch,  Dunbar  Laboratory 

Dr.  George  P.  Carrier,  309  Pierce  Hall 

Illinois  Institute  of  Technology 
Chicago,  Illinois 

Attn:  Dr.  L.  H.  Donnell 

Dr.  E.  Sternberg 
Dr.  W.  Osgood 

University  of  Illinois 
Urbana,  Illinois 

Attn:  Dr.  N.  M.  Newmark 

Ecif;ineering 
T.  J.  Dolan 

Dr.  F.  Seits,  Department  of  Physics 
Departmont  or  Theoretical  and  Applied  Mathema- 
tics, Attn:  Prof.  M.  C.  Steale 

Indiana  University 
Bloomin''ton,  Indiana 

Attn:  Dr<  T.  Y.  Thomas 

Institute  for  Advanced  Study 
Princeton,  New  Jei'sey 

Attn;  Dr.  J.  vcn  Nexunann 

Iowa  State  College 
Ames.  Iowa 

Attn:  Dr.  G.  Murphy 

Dr.  D.  L.  Roll 


Johns  Uopkins  University 
Baltimore,  Maryland 

Attn:  Dr.  W.  H.  Hoppman,  II 

Director,  Applied  Physics  Laboratory 
Johns  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 


Lehigh  University 
Bethlehem,  Pennsylvania 

Attn:  Mr.  Lynn  S.  Beedle 


Jfessachusetts 
Cambridge  39, 
Attn:  Dr. 

Dr. 
Dr. 
Dr. 
Dr. 
Dr. 
Dr. 
Dr. 
Dr. 


Inr,titute  of  Technology 

I^assachusetts 

F.  B.  Hildebrand 

J.  M.  Lessels 

V7.  M.  Ifurray 

E.  Relssner 

H.  Cohen,  Rm.  8-413,  Dept,  of  Metallurgy 
B.  L.  Averbach,  Department  of  Metallurgy 
J.  T.  Norton 
E.  Orowan 

R.  Bisplinghoff , Dept.  Aero.  Engr. 


M-1 

M-1 

R-1 

R-1 


R-1 

H-1 


R-1 


R-1 


R-1 


R-1 

R-1 

R-1 


R-1 

R-1 

R-1 


M-1 


R-1 

R-1 

R-1 

R-1 


M-1 


6 


University  of  Michigan 
Ann  Arbor)  Michigan 

Attn:  Dr*  Bruce  G«  Johnston 

Dr*  Paul  Nagdhl)  Dept*  of  Engineering  Mechanics 
Dr*  N*  Coburn 
Dr*  W*  Kaplan 

New  York  University 

Institute  for  Mathematics  & Mechanics 

45  Fourth  Avenue 

New  York  3,  Now  York 

Attn:  Professor  R.  Courant 

Dr*  G*  Hudson 

New  York  University 
Now  York  53 > New  York 

Attn:  Dr.  C.  T.  Wang)  Dept*  of  Aeronautics 

Northwestern  University 
Evanston,  Illinois 

Attn:  Dr*  M*  M*  Hetenyl 

lAilverslty  of  Notre  Daise 
Notre  DamO)  Indiana 
Attn:  Dr*  P.  A.  Beck 

Pennsylvania  State  College 
State  College)  Pennsylvania 
Attn:  Dr*  J*  A*  Sauer 

Dr*  Joseph  Marin 
Dr*  J*  W*  Fredrickson 

Princeton  University 
Princeton,  New  Jersey 
Attn:  Dr*  S*  Lefschets 

Dr*  L*  Lees 
Dr*  J*  V*  Charyk 

Purdue  University 
Lafayette,  Indiana 

Attn:  Dr*  C*  A*  Erlngen 

Rensselaer  Polytechnic  Institute 
Troy.  New  York 
Attn:  Library 

Dr*  Paul  Leiber 
Dr*  C*  0*  Dohrenwond 
Dr*  G*  H*  Lee 

Santa  Clara  University 
Santa  Clara,  California 
Attn:  Dr*  R*  M,  Hermes 


7 


R-X 

H-X 

R-X 


R-X 


M-X 

R-X 


R-X 


M-X 


R-X 


R-X 


R-X 


M-X 

M-X 


R-X 


R-X 

R-X 

R-X 


R-X 


Stanford  Univarslty 
Stanford,  CaXlfornla 
Attn;  Dr*  L*  Jacobsen 

Dr*  A*  PhlXXlps,  Dept*  of  Mechanic aX  Engineering 
Dr*  J*  N*  Goodler 

Stevens  Institute  of  Technology 
Hoboken,  New  Jersey 

Attn;  Dr*  E*  G*  Schneider 

Swarthmore  College 
Swarthmore,  Pennsylvania 
Attn:  Capt*  W.  P*  Roop 

Dr*  S.  T.  Carpenter 

University  of  Texas 
Austin  12,  Texas 

Attn;  Dr*  A*  A*  Topractsoglou 

University  of  Utah 
Salt  Lake  City,  Utah 
Attn:  Dr*  H*  ^rlng 

Washington  State  College 
Pullman,  Washington 
Attn;  Dr*  B*  Fried 

Wheaton  College 
Norton,  Massachusetts 
Attn:  Dr*  H*  Gelrlnger 

Aerojet,  Inc* 

Azusa.  California 
Attn;  P*  Zwicky 

Aluminum  Company  of  America 
New  Kensington,  Pennsylvania 
Attn;  R*  L*  Templin 

H*  N*  Hill,  Aluminum  Research  Laboratory 

Armstrong  Cork  Company 
Lancaster,  Pennsylvania 
Attn;  J«  W*  Scott 

Bell  Telephone  Laboratories 
Murray  Hlll^  New  Jersey 
Attn:  C*  Herring 

D.  P*  Ling 
W*  P*  Mason 

Corning  Glass  Company 
Corning,  New  York 

Attn;  J*  T*  Littleton 


8 


E.  I*  Dupont  de  Nemoxirs  & Co* , Inc* 

Wilmington  98,  Delaware  ^ ^ e 

R-1  Attn*  J.  H.  Paupol,  Materials  of  Construction  Section 

General  Electric  Company 
Schenectady,  New  York 
R-1  Attn:  H*  Fehr 

R-1  H«  Poritsky 

R-1  J.  H*  Hollomon 

General  Motors 
Detroit,  Michigan 
R-1  Attn:  J.  0.  Almen 

Lockheed  Aircraft  Company 
Department  72-25,  Factory  A-1,  Building  66 
Burbank,  California 
R-1  Attn*  Engineering  Library 

Midwest  Research  Institute 
Kansas  City,  Missouri 
R-1  Attn:  M.  Goland 

Pratt  & IVhltney  Aircraft  Corporation 
East  Hartford,  Connecticut 
R-1  Attn:  R.  Morrison 

U.  S.  Rubber  Company 
Passaic,  New  Jersey 
R-1  Attn:  H.  Smallwood 

Welding  Research  Council 
Engineering  Foundation 
29  V/est  39  Street 
New  York  18,  New  York 
M-1  Attn:  W#  Spraragen,  Director 

Westlnghouse  Research  Laboratories 
East  Pittsburgh,  Pennsylvania 
R-1  Attn:  Dr*  £.  A*  Davis 

M-1  Dr.  A.  Nadal 

136  Cherry  Valley  Road 
Pittsburgh  21,  Pennsylvania 

Westlnghouse  Electric  Corporation 
Lester  Branch  P«0« 

Philadelphia,  Pennsylvania 

R-1  Attn:  R.  P.  Kroon,  Manager  of  Englneering,AGT  Division 

University  of  Pennsylvania 
Towne  Scientific  School 
Philadelphia  ^•,  Pennsylvania 

R-1  Attn*  Dr.  C.  W.  flacGregor,  Vice  President  in  Charge  of 

Scientific  and  Engineering  Studies 


R-1  Professor  G.  Wastlund 

Cement  & Concrete  Research  Institute 
Royal  Institute  of  Technology 
Stockholm  70 f SWEDEN 

M-1  Professor  John  E.  Goldberg 

Department  of  Structural  Engineering 
Purdue  University 
Lafayette,  Indiana 

R-1  Jttlad  F.  Hanna 

Jg  P.S.S.B. 

Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts 

R-1  Dr.  W,  Frelberger 

Department  of  Supply 
Aeronautical  Research  Laboratories 
Box  4331  GPO 
Melbourne,  AUSTRALIA 

M-1  Professor  B.  W.  Shaffer 

Department  of  Mechanical  Engineering 
New  York  University 
New  York  53,  New  York 

M-1  Professor  G.  Sachs 

Division  of  Metallurgical  Research 
Engineering  & Science  Campus 
East  Syracuse  4,  N,  Y, 

R-1  Professor  Asia  Ghali 

Head.  Structural  Analysis  Department 
Fonaa  University 
Giza,  EGYPT 

M-1  Professor  Hugh  Ford 

Mechanical  Engineering  Department 
Imperial  College  of  Science  & Technology 
London,  S.W.7 
ENGLAND 

R-1  Dr.  R.  H.  Wood 

D.  S.  I.R«  Building  Research  Station 

Garston,  Watford,  Herts 

ENGLAND 


